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How many supernovae are we missing at high redshift? 
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ABSTRACT 

Near-infrared and radio searches for core-collapse supernovae (CC SNe) in the local 
universe have shown that the vast majority of the events occurring in massive starburst 
are missed by the current optical searches as they explode in very dusty environments. 
Recent infrared observations have shown that the fraction of star-formation activity 
that takes place in very luminous dusty starbursts sharply increases with redshift and 
becomes the dominant star formation component at z^0.5. As a consequence, an in- 
creasing fraction of SNe are expected to be missed by high-redshift optical searches. 
We estimate that 5-10% of the local CC SNe are out of reach of the optical searches. 
The fraction of missing events rises sharply toward z=l, when about 30% of the CC 
SNe will be undetected. At z=2 the missing fraction will be about 60%. Correspond- 
ingly, for type la SNe, our computations provide missing fractions of 15% at z=l and 
35% at z=2. Such large corrections are crucially important to compare the observed 
SN rate with the expectations from the evolution of the cosmic star formation history, 
and to design the future SN searches at high redshifts. 
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1 INTRODUCTION 
1.1 Infrared galaxies 

Most of the star formation activity in the local universe oc- 
curs in "normal" galaxies characterized by star formation 
rates (SFRs) of 0.1-5 M0/yr and moderate dust extinctions 
(Av < 1). Nevertheless, a few objects with high SFRs (100- 
1000 Mo/yr) and heavy dust attenuation (Ay =5-50 mag) 
exists. Only a very small fraction of the UV light produced 
by their young stars can escape these galaxies, while most 
of the energy is converted into thermal infrared (IR) radi- 
ation by heating the dust (e.g., Burgarella et al. 2006; see 
also Mannucci & Beckwith, 1995). As a consequence their 
bolometric luminosity is dominated by the far-IR radiation 
from the hot dust. These galaxies are classified as Luminous 
IR Galaxies (LIRGs) if their total far-IR luminosity is above 

10 11 L Q (corresponding to SFR~10-20 M /yr, Choi et al. 
2006) and as Ultraluminous IR galaxies (ULIRGs) above 

10 12 L0 . Most LIRGs and ULIRGs show the signs of recent 
major mergers (e.g., Lagache et al. 2005) that are consid- 
ered to be the origin of the strong episode of star formation. 
Liang et al. (2004), Marcillac et al. (2006) and Choi et al. 
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(2006) have studied samples of distant (z~0.8) LIRGs and 
found no evolution of the properties of these galaxies with 
redshift. 



Even if these objects are very active in forming stars, 
their comoving volume density is low in the local universe 
and their contribution to the total star formation activity 
is below 5% (Soifer & Neugebauer, 1991). Two decades 
ago, using data from the IRAS satellite, it was discovered 
that massive starburst galaxies were much more numerous 
and/or luminous in the past than they are today. For 
example, Hacking et al. (1987) have shown that the counts 
of the bright IRAS galaxies cannot be explained without a 
fast evolution up to z=0.2. Using the Spitzer satellite, the 
evidence for rapid evolution has been recently extended up 
to z=l (Le Floch et al. 2005) and beyond (Perez-Gonzalez 
et al. 2005; Daddi et al. 2005). As a consequence, it is 
now well established that massive dusty starburst are rare 
objects in the local universe but were much more common 
a few Gyr ago and dominated the star formation (SF) 
activity before z=2. 



2 Mannucci, Delia Valle & Panagia 



1.2 Core-Collapse supernovae 

In the last few years, another class of astronomical objects, 
supernovae, (SNe) has proved to be very useful in observa- 
tional cosmology. Besides the well known results about the 
accelerated expansion of the universe obtained from studies 
of type la SNe, also CC SNe in high-redshift galaxies have at- 
tracted much interest. These objects, which are generated by 
the gravitational collapse of very massive stars (M>8M0), 
are among the most important producers of metals (e.g., 
Matteucci & Greggio, 1986), are considered to be responsi- 
ble for dust production in the young universe (Maiolino et 
al. 2004a, 2004b), could dominate the feedback in the pro- 
cess of galaxy formation (see, for example, Ferrara & Ricotti, 
2006, and references therein), and can be used to estimate 
the SF density at high redshifts (e.g., Dahlen et al. 2004). 

The measured rates of CC SNe both at low and at high 
redshifts are based on optical observations because only at 
these wavelengths the current instrumentation has sufficient 
field coverage, spatial resolution and sensitivity to detect 
large numbers of SNe within a reasonable observing time. 
In the local universe (z<0.1), the best available rates have 
been computed by Cappellaro et al. (1999) and Mannucci et 
al. (2005) and are based on a compilation of a few visual and 
photographic searches. Other searches for CC SNe based on 
CCD imaging are available or currently active (Filippenko 
et al. 2001; Strolger, 2003). Recently, the rates for higher 
redshifts (z~0.5) CC SNe have been measured by Dahlen 
et al. (2004) and Cappellaro et al. (2005), although large 
uncertainties remain, for example because of the lack of an 
extensive spectral identification of the SNe candidates. To 
derive the intrinsic SN rates from the observed number of 
events it is necessary to estimate the fraction of SNe that 
are missed because of any known reason, as dust extinction 
in the host galaxies and the reduction of detection efficiency 
near the galaxy nucleus. Customarily, these corrections are 
estimated in the local universe and are assumed to hold un- 
changed with redshift. Also, they are appropriate to correct 
the SN rates in normal galaxies and do not take into ac- 
count the presence of LIRG and ULIRG among the galaxy 
target. For example, Dahlen et al. (2004) assume a mod- 
erate average extinction of E(B-V)=0.15, corresponding to 
Ay ~ 0.5, appropriate for the local population of normal 
galaxies, and use this correction at all redshifts. The result- 
ing rate accounts only for the SNe exploding in galaxies 
where the star formation activity occurs in relatively clean 
enough environments so that significant fraction of SNe can 
be detected. SNe exploding inside dusty starbursts are miss- 
ing from these counts. 



1.3 Supernova searches at longer wavelengths 

In the local universe, optical searches for SNe in LIRGs and 
other starburst galaxies (Richmond et al. 1998) failed to de- 
tect the expected enhancement of the rate over more quies- 
cent galaxies due to the strong activity of star formation. To 
test the hypothesis that this shortage is due to the presence 
of large amounts of dust obscuring a substantial fraction of 
SNe, several groups (van Buren & Norman, 1989; Maiolino 
et al. 2002; Mannucci et al. 2003; Mattila et al. 2004; Cresci 
et al. 2006) have carried out near-IR searches for SNe in 
starburst galaxies. As dust extinction is about 10 times less 



important in the K-band than in the optical, these searches 
were expected to detected substantially higher rates. Indeed, 
some of these searches have detected a population of near-IR 
CC SNe that are missed by optical searches. For example, 2 
SNe (SN1999gw, Cresci et al. 2002, and SN2001db, Maiolino 
et al. 2001) out of the 4 events detected in the K-band by 
Maiolino et al. (2002) are missing from the optical catalogs. 
At least one of them is characterized by 5-6 magnitudes of 
extinction in the V band which would have prevented its 
detection at optical wavelengths. However, these searches 
failed to recover the majority of the expected events. Man- 
nucci et al. (2003) demonstrated that about 80% of the CC 
SNe expected to explode in LIRG and ULIRG, on the basis 
of their far-IR luminosity, are still undetected even by near- 
IR searches. Similar results have been obtained by Cresci et 
al. (2006) for a SN search with the near-IR camera NICMOS 
on the HST. The shortage of SNe can be due either to the 
presence of very large dust column densities (Ay > 20) or 
to the fact that most of the SNe explode in the nucleus of 
the starburst galaxy, where SNe cannot be easily disentan- 
gled from the bright emission of the host galaxy. In both 
cases, most of these SNe are missing from the estimates of 
the local rates, and their rate cannot be inferred from the 
SNe detected in normal galaxies. Recently, Lonsdale et al. 
(2006) have monitored the nucleus of Arp 220, the prototypi- 
cal ULIRG, at radio wavelengths with sub-arcsec resolution. 
They detected a few faint variable sources and a population 
of brighter, slowly variable point sources. If these data are 
interpreted in terms of SNe and SN remnants, the number 
of SN events would be enough to power all the far-IR lu- 
minosity of Arp 220. However, none of this object could be 
detected by optical or near-IR searches as they are deeply 
embedded in the huge amount of nuclear dust. 

The effect of dust extinction over the observed SN 
rate is expected to be severe for LIRG and ULIRG. Many 
authors, such as Cram et al. (1998), Jansen et al. (2001), 
Hopkins et al. (2001), Sullivan et al. (2001), Kewley et 
al. (2002), Perez-Gonzalez et al. (2003) and Choi et al. 
(2006), have studied the amount of extinction in starburst 
galaxies as a function of the SFR, and concluded that 
the two quantities are strongly correlated, i.e., the more 
active galaxies also show higher levels of dust extinction. 
As a consequence, both the SFR and the SN rate can be 
derived by optical and UV observations for galaxies with 
SFR^5M0/yr because the low average extinction (Av ^ 1) 
allows for a meaningful correction. On the contrary, for 
larger SFRs, the average extinction becomes very large. 
For example Choi et al. (2006) have studied a sample of 
LIRG at z^0.8 and have obtained an average extinction of 
Av = 2.5, which corresponds to an average reduction of the 
flux in the [OII]A3727 line of 97%. 

For these galaxies the total SF activity can only be derived 
at wavelengths less sensitive to dust extinction (far-IR and 
radio), while the total SN rate cannot be measured. 

In this paper we want to estimate the fraction of SNe 
which are missed both in the local and in the distant uni- 
verse by current and future optical or near-IR SN searches. 
In the next section we will review the information on the 
SF activity at z<2 as derived from UV and far-IR obser- 
vations in order to quantify the star-formation activity oc- 
curring in clean and dusty environments. In section [3] we 
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Figure 1. Left panel: Evolution of the SF density with redshift as derived from UV observations, uncorrected for dust extinction. Data 
points are from Lilly et al. (1996), Connolly et al. (1997), Treyer et al. (1998), Cowic ct al. (1999), Sullivan et al. (2000), Massarotti ct 
al. (2001), Wilson et al. (2002), Wolf et al. (2003), Schiminovich et al. (2005), Baldry et al. (2005), Wyder et al. (2005), Thompson et al. 
(2006), and Dahlen et al. (2006). Original data have been first converted into luminosity density at 2800A using a UV slope (3 = —1.5, 
and then to a SF density as in Baldry et al. (2005). The solid line is the fit to these data points in eq. 2] Right panel: Evolution of the SF 
density with redshift as derived from IR and radio observations. The data points are from Condon (1989), Flores et al. (1999), Barger 
et al. (2000), Haarsma et al. (2000), Machalski & Godlowski (2000), Sadler et al. (2002), Serjeant et al. (2002), Condon et al. (2002), 
Georgakakis et al. (2003), Perez-Gonzalez et al. (2005), Mauch (2005) and Babbedge et al. (2006). The dashed regions in green, orange 
and red show the range of possible contributions from FIRG, LIRG and ULIRG, respectively, as derived by Perez-Gonzalez et al. (2005). 
The thick lines in the same colors are the assumed evolution of the three populations as described by eqs.[T]|2] and [3] The magenta line 
is the total SF density given by the sum of the three contributions, and can be directly compared with the data points. 



will use the resulting SF density to derive the SN rate. 
Throughout this paper we will use the concordance cosmol- 
ogy (/iioo,n m ,n A ) = (0.7,0.3,0.7). 



2 EVOLUTION OF THE LUMINOUS 
INFRARED GALAXIES 

In the last decade, many authors have contributed to plot 
a coherent picture of the evolution of the star formation 
activity along the Hubble time (see Hopkins & Beacom 
2006 for a recent review). Observations at radio, IR, optical, 
UV and X-ray wavelengths have detected the SF activity 
over a broad range of redshifts. The emerging picture shows 
a fast increase of the SF density with increasing redshift 
between and 1, a roughly constant density between 
z~l and z~4, and a possible decrease for z>4 (Bouwens 
& Illingworth 2006, Mannucci et al. 2007). Observations 
made at different wavelengths sample different aspects 
of the star formation activity. The UV observations are 
deeply affected by the presence of dust, and evidence 
of significant extinction is found in practically all types 
of high-redshift star-forming galaxies. On the contrary, 



far-IR observations are more suitable to trace the star 
formation that takes place in dusty environments. It is 
now clear that a complete census of the star formation can 
be obtained only by putting together all these contributions. 

Perez- Gonzalez et al. (2005) have used Spitzer data at 
24/^m of about 8000 galaxies with reliable spectroscopic or 
photometric redshifts to constraint the evolution of the far- 
IR Luminosity function (LF) with redshift. They confirmed 
the strong evolution of the far-IR luminosity density with 
redshift and, most importantly, were able to trace the dif- 
ferent contributions to this evolution (see Fig. 1). At z<l, 
the dominant contribution comes from galaxies with low IR 
luminosities (faint IR galaxies, FIRGs, L< IO^Lq). The 
LIRG component accounts for a few percent of the total 
SF density in the local universe but has a faster evolution 
with redshift, becoming the dominant component at z~l. 
At higher redshifts, the LIRG contribution declines and the 
SF density is dominated by ULIRGs. 

Several effects contribute to the uncertainties of this pic- 
ture, especially at z>l: the value of the faint-end slope a of 
the LF is poorly constrained; the total far-IR luminosity in 
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the 8-1000/mi wavelength range is estimated from the flux 
at one single wavelength; most of the redshift determinations 
are photometric, allowing the presence of a significant frac- 
tion of interlopers. Despite all these uncertainties, the fast 
evolution up to z—1 is a strong results confirmed by many 
studies, as explained in the introduction. For example, us- 
ing a smaller sample of galaxies but with a larger fraction 
of spectroscopic redshifts, Le Floch et al. (2005) found the 
same type of evolution for FIRG, LIRG and ULIRG up to 
z=l, even if they found a larger total number of LIRG. Also 
Chary et al. (2007) find the same behavior, with a "FIRG 
epoch" at z<l, a LIRG epoch at l<z<2, and a ULIRG epoch 
at z>2. The most uncertain effect is the reduction in num- 
ber of LIRG at z>l and the corresponding increase in the 
ULIRG comoving density, but this behavior has hardly any 
effect on the results of this paper because the fractions of 
SNe detected in LIRG and ULIRG are very similar (see next 
section) . 

Fig. 1 shows the observed evolution of the SF density 
with redshift as measured in the UV and in the far-IR. Data 
points are derived from the sources listed in the figure cap- 
tion and assume the "SalA" modified Salpeter initial mass 
function (IMF, Baldry & Glazebrook 2003, see also Hopkins 
& Beacom 2006). UV data are uncorrected for dust extinc- 
tion. The dashed areas in the IR panel show the uncertain- 
ties of the various contributions from Perez-Gonzalez et al. 
(2005). The solid lines show the SFD deriving from the ob- 
served UV galaxies, FIRG, LIRG and ULIRG. To describe 
these curves in an analytic form, we use a parametrization in 
which a power-law of (1+z) is multiplied by an exponential 
cutoff at high redshifts: 



\og{SF F1RG ) 



log(SF LIRG ) 



log(SF ULIB .G) 



log(SF v 



= -2.05 + 3.5 log(l + «) 

-log [l + cxp(3.7 (z- 1.1)] 

= -3.65 + 8.8 log(l + z) 

- log [1 + exp(4.6 (z - 0.96)] 

= -5.90 + 12.5 log(l + z) 

- log [1 + cxp(3.0 (2-1.1)] 

-2.20 + 2.6 log(l + z) if z< 1 
-1.43 if z > 1 



(1) 

(2) 

(3) 
(4) 



These analytic equations are valid up to z=2. Since 
these contributions to the cosmic SF activity occurs in dif- 
ferent environments, each of them is associated to a differ- 
ent fraction of SNe that can be detected. As a consequence, 
while the evolution of the intrinsic CC SN rate follows the 
total SF density, the observed SN rate is expected to show 
a different behavior, because the fraction of detectable SN 
changes with redshift. 



3 INTRINSIC AND OBSERVED SN RATES 
3.1 Core-Collapse Supernovae 

The total CC SN rate at any redshift is proportional to the 
total SFD at that epoch as the lifetime of the exploding 
stars is short, of the order of 10 7 yr, and can be neglected. 
The proportionality between SN rate and SFD depends on 



the assumed IMF. If p is the SN rate in units of SN/yr and 
ip is the SFR in units of Mo/yr, a simple integration of the 
Salpeter IMF between 8 and 50 M Q provides p = 0.007i/> 
(see, for example, Cappellaro et al. 2005). For the SalA 
IMF, Hopkins & Beacomp (2006) obtain p = 0.00915^. 
Nevertheless, it should be noted that the total SN rate 
expected from the UV and far-IR observations is practically 
independent of the assumed IMF, as both the UV and 
far-IR luminosities are mostly due to the same massive 
stars producing the CC SNe. For example, the use of a 
Salpeter instead of a SalA IMF would result in a SFD 
30% higher than that shown in Fig. 1, but the expected 
number of SNe for a given SFR would be 30% lower, 
and the two factors would cancel out each other. On the 
contrary, the stellar mass range to produce a CC SNe is 
an important parameter which is only weakly constrained 
by models. Here we use the "canonical" range between 8 
and 50 MQ(e.g., Woosley & Weaver, 1986), while the use of 
different mass ranges would produce different total number 
of CC SNe. For example, using a range between 10 and 40 
M would reduce the expected rate of 32%. The total CC 
SN rate we derive is shown in Fig. 2. It is very similar to 
that in Hopkins & Beacom (2006) apart from the different 
parametrization used for the SFD. 

To derive the observed SN rate from the measured SF 
density we need to add a few assumptions about the average 
properties of the galaxies of the different populations: 

(i) Since the star formation activity detected in UV se- 
lected galaxies and in FIRG occurs in relatively clean envi- 
ronments, we assume that the optical SN searches are ade- 
quate to measure the rates. This means that many CC SNe 
in these galaxies are detected, and that the missing part can 
be recovered by estimating the average dust extinction as 
in Cappellaro et al. (1999) or in Dahlen et al. (2004). This 
is justified by the fact that low values of Ay are usually 
found in these galaxies. Also at high- redshift, most of the 
UV-selected Lyman Break Galaxies (LBG) show Ay < 1 
(Vijh et al. 2003; Burgarella et al. 2006). For these reason, 
we assumed that all the SNe in UV-selected galaxies and at 
least 80% of those in FIRG are accounted for by the optical 
searches, either directly or through the extinction correc- 
tions. 

(ii) In the local universe, less that 10% of the SNe in 
LIRG are detected by optical searches, and less then 20% 
can be recovered by near-IR searches (Mannucci et al. 2003). 
At higher redshifts, the fraction of detected SNe could be 
higher if distant LIRG have a much smaller content of dust, 
for example due to a lower metallicity. To allow for any pos- 
sible, although unobserved, evolution in the dust properties 
of LIRG, we assume that the maximum fraction of SNe that 
are accounted for by the optical searches is 0.2 at z=0 and 
linearly increases to 0.4 at z=2. This is a solid upper limit 
as optical searches are not likely to detect more that 40% of 
the SNe in LIRG at any redshift. The minimum fraction of 
detectable SNe is assumed to be at any redshift. 

(iii) Only a minor fraction of SNe are detected in the 
ULIRG in the local universe, as these galaxies are known to 
enshroud their star formation activity inside large amounts 
of dust. For example, Arp 220 is expected to produce 4 SNe 
per year (Mannucci et al. 2003), while only one possible SN 
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Figure 2. Le/4 panel: Evolution of the CC SN rate with redshift. The dashed line shows the intrinsic total SN rate derived from the 
SFD. The SN rate that can be recovered by optical searches is shown in light grey, and is compared with the data points from Cappellaro 
et al. (2005) (squares), and Dahlen et al. (2004) (circles). For consistency, the former two data points have been converted to comoving 
rates as in the latter work. In dark grey we show the expected SN rate for SNe that cannot be detected because they are exploding in 
very dusty star burst, as explained in the text. The reported uncertainties correspond to the possible range of the fraction of SNe that 
can be detected in the various classes of galaxies, as explained in Sect. [3] Right panel: Fraction of CC SNe than are not present in the 
optical and near-IR searches as a function of redshift. The dashed line shows the results of eq. [5] 



was ever detected (Cresci et al. 2006), and only by using the 
HST near-IR camera. As a consequence, we assume that 
between and 10% of the SNe in ULIRG are present in the 
rates derived from optical/near-IR observations. 



All these quantities are rather uncertain, but the assumed 
minimum and maximum fractions are expected to encom- 
pass all the possible values. 

By using these assumptions and the SFDs described in 
eqs.[T}{4]we derive the fraction of SNe that can be recovered 
by optical searches. The results are shown in Fig. 2 where 
the expectations are compared with the current set of ob- 
servations. In the local universe only a small fraction of SNe 
(5-10%) are not accounted for by the optical rates, a frac- 
tion which is smaller that the uncertainties in the measured 
rate. Going toward higher redshifts, this fraction increases 
rapidly, reflecting the rapid evolution of the number density 
of LIRGs and ULIRGs. At z=l the undetected fraction is 
20-40%, rising to ~60% at z=2. 

The behavior of the missing fraction /(CC) of CC SNe 
can be described by a linear relation with redshift: 



/(CC) = 0.05 + 0.28 



for z < 2 



(•») 



3.2 Type la supernovae 

Type la SNe are also affected by dust extinction as they 
explode in all environments. For example, in the local uni- 
verse there are type la SNe that show large amounts of 
dust extinction, such as SN2002cv (Di Paola et al. 2002) 
and SN2000E (Valentini et al. 2003). Chary et al. (2005) 
have studied 50 la SN host galaxies at 0.1<z<1.7 and found 
larger amounts of dust than in normal field galaxies, and far- 
IR luminosities similar to those in CC SN hosts. This is an 
indication that the environments of the two SN populations 
are similar. 

While all the CC SNe are expected to explode soon after 
the formation of the massive parent star, the la SNe show 
a wide delay time distribution (DTD; i.e., a wide distribu- 
tion of the time elapsed between the birth of the progen- 
itor star and the SN explosion). According to many com- 
putations (e.g., Greggio & Renzini 1983, Greggio 2005) the 
DTD could range between 10 7 and 10 10 years. Such a large 
spread is needed to reproduce the variation of the la SN rate 
per unit mass with host galaxy color (Mannucci et al. 2005, 
2006; Greggio, 2005) and to account for the metallicity evo- 
lution of the galaxies (e.g., Matteucci et al. 2006, Calura & 
Matteucci 2006). Mannucci et al. (2006) have shown that 
the dependence of the la SN rate in early-type galaxies with 
the radio power of the host (Delia Valle et al. 2005) is best 
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Figure 3. Left panel: Evolution of the type la SN rate with redshift. See the caption of Fig. 2 for the explanation of the symbols. The 
observed rates, shown by the dots, are from Mannucci et al. (2005), Magdwick et al. (2003), Strolger (2003), Blanc et al. (2004), Pain 
et al. (2002), Neill et al. (2006) and Dahlen et al. (2004). Right panel: Fraction of type la SNe than are not present in the optical and 
near-IR searches as a function of redshift. The dashed line shows the results of eq. [6] 



explained if about half of the la SNe, the so called prompt 
component, explode within 10 8 yr from star formation, while 
the rest (the tardy component) has a much wider DTD. For 
this reason about half of the la SNe (the prompt component) 
should have extinction properties similar to the CC SNe as 
their progenitors explode on similar timescales. The time 
evolution of the extinction of the tardy component is more 
uncertain. In a naive picture we can assume that all the 
prompt population suffers of the same extinction of the CC 
SNe, as explained in section 13.11 while the tardy population 
has time to escape out of the very dusty star-forming region 
and, as a consequence, show a much milder extinction and 
can easily be recovered by observations. Figure 3 shows the 
results of the computations based on these hypothesis and 
using Mannucci et al. (2006) DTD: the fraction of missing 
la SNe increases from a few percent in the local universe to 
about 15% at z=l and 35% at z=2. This behavior can be 
described by: 



/(la) =0.02 + 0.12 z 1 



for z < 2 



(6) 



We point out that the results for type la SNe are somewhat 
more uncertain than those for CC SNe, because they also 
depend on the actual DTD and on the evolution of the 
extinction with delay time. 

The expected volumetric rates in Figures 2 and 3 can be 
translated into the expected number of SNe per sq. arcmin 
per year that can be detected by any survey as a function 
of the maximum redshift reached, as plotted in figure 4. For 



example, a survey capable of detecting CC SNe up to z=2 
is expected to find 3 events per year for each sq. arcmin 
monitored out of the 5 CC SNe exploded. 



4 DISCUSSION 
4.1 Uncertainties 

The results in Fig. 2 and Fig. 3 depend on the various as- 
sumptions, especially the somewhat arbitrary division of the 
galaxies in FIRG, LIRG and ULIRG, and the fraction of SNe 
that are assumed to be detected in each of these classes. 
Even if some uncertainties are present, a robust conclusion 
is inevitable: as most of the the star formation activity at 
z>l occurs in dust- rich LIRGs and ULIRGs, optical SN 
searches will not follow the global SN rate but only the 
fraction of detectable SNe, which will decrease sharply as 
the redshift increases. In particular, it is impossible that the 
fraction of SNe missed at high redshift is similar to that lost 
in the local universe (about 10% for the CC SNe). For ex- 
ample, such a low level of incompleteness would be possible 
at z=l only if more than 75% of the CC SNe in LIRGs are 
dected by optical/near-IR searches; and at z=1.5, not only 
all the CC SNe in LIRG should be detected, but also 40% 
of those in ULIRGs. These high detection efficiencies would 
imply a strong evolution of the dust properties of LIRGs and 
ULIRGs, evolution that is not observed (Liang et al. 2004, 
Marcillac et al. 2006, Choi et al. 2006). As a consequence, 
although the exact amount remains uncertain, we conclude 
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Figure 4. Predicted cumulative number of CC and la SNe ex- 
ploding in one year in 1 sq. arcmin of sky from z=0 up to a given 
redshift. The dashed line shows the total number of exploding 
SNe, the solid line the SNe that can be detected by optical and 
near-IR, searches. 

that at z>l a large fraction of SNe remain unaccessible to 
optical/near-IR searches. 

4.2 Comparison with the observed rates for 
Core-Collapse SNe 

Hopkins & Beacom (2006) have already derived the total CC 
SN rate from the cosmic SFD. In their Fig. 7 they compared 
the expectations with the observations, and found that in all 
cases the observed rates are below the expectations. They 
attribute this finding to several effects, such as incomplete- 
ness of the local rate and insufficient dust correction for 
distant rates. The same effect can be seen in our Fig. 2, 
but in this case the comparison must be done with the ob- 
served fraction instead of the total one. The discrepancy is 
much reduced, especially at high redshifts: the presence of a 
significant fraction of SNe that cannot be detected helps in 
matching the expected and observed SN rates. 

This remaining systematic overprediction of the rates 
can be attributed to several effects. As discussed in section|3l 
it could be due the mass range for CC production, which 
could be significantly different (narrower) than the 8-50 Mq 
range as adopted here. 

Another potentially important effect is related to the 
origin of the far-IR flux in the FIRG galaxies. It is well 
known that in the quiescent galaxies a significant fraction of 
far-IR luminosity is not linked to the young stellar popula- 
tions but rather to the general radiation field. For this reason 
the SFD for the FIRG in figure 1, derived from the far-IR 
data, could be overestimated because part of this far-IR flux 
could be related to old rather than young stars. As the SN 
counts in the local universe are dominated by this FIRG 
component, the actual total number of SN at low redshift 



could be lower than the predictions shown in the left panel 
of figure 2, more in agreement with the measured rates. At 
variance, the fraction of missing SN in the right panel of fig- 
ure 2 is not very sensitive to the total normalization of SN 
in FIRG as it effects both detected and undetected SNe. 

It should be noted that the local rate from Cappellaro 
et al. (1999) (0.59 ± 0.22 ■ 10" 4 Mpc" 3 yr _1 ) is the only 
value more then la below the expectations (1.3 ± 0.1 ■ 10 -4 
Mpc -3 yr _1 ). Indeed, this value is based on a compilation 
of visual and photographic searches that could miss a 
significant fraction of SNe, for example near galaxy nuclei. 
Cappellaro et al. (1999) considered this effect, but the 
applied correction may have been not large enough. This 
fact can be understood considering the recent results of 
the Lick Observatory SN Search (LOSS, Filippenko et al., 
2001). The LOSS search has a deeper limiting magnitude 
and, in practice, is a volume limited search also for the 
CC SNe. The fraction of CC SNe in the LOSS sample 
(Li, 2005) is considerably larger than in Cappellaro et al. 
(1999) sample, indicating that the latter sample may miss 
a significant fraction of CC SNe. 



4.3 Comparison with the observed rates for type 
la SNe 

Many more measurements of type la SN rate are available 
(see Fig. 3) although the spread among different results is 
sometimes larger than the quoted errors, a fact that sug- 
gests the presence of uncorrected systematic effects. Data 
show a marked increase of the rate between z=0 and z=l 
and, possibly, a constant rate at higher redshifts. The wide 
DTD reproducing the rates in the local universe in the var- 
ious classes of galaxies (for a full discussion see Mannucci 
et al., 2006) tends to produce more SNe at z>l than ob- 
served by Dahlen et al. (2004). Even if these observations 
are very uncertain, as they are based on a small number 
of confirmed SNe and on several untested assumptions, the 
presence of dust extinction helps in reconciling expectations 
with observations. 

4.4 Other sources of incompleteness 

In this paper we have only discussed the incompleteness of 
the measured rates due to the increasing fraction of SFD 
occurring in dusty environments. There are other possible 
sources of incompleteness that could be present and could 
give important contributions at high redshift. For example, 
a class of faint CC SNe have been recently detected in the 
local universe (Pastorello et al. 2004) that could be associ- 
ated with explosion energy so small that most of the 56 Ni 
falls back onto the compact stellar remnant. Such an ef- 
fect was also invoked by Delia Valle et al. (2006) to explain 
the absence of any observable SN associated with the long 
GRB060614. These massive stars contribute to the far-IR 
luminosity of the galaxies but not the observed CC SN rate. 
The objects of this class seem to be quite rare in the local 
universe, about 5% of the total type II rate according to 
Pastorello et al. (2004) , but the evolution of their rate with 
with redshift is not known. If they become more common 
at high redshifts, for example due to metallicity effects, the 
expected rates will be lower than presented here. 
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Summarizing, on the basis of the recent estimates of 
the density evolution of the star formation activity in clean 
and dusty environments, we have estimated the fraction of 
SNe that are likely to be missed by the optical searches up 
to z=2. We find a strong evolution of the missing fraction: 
for CC SNe it is as small as 5-10% in the local universe and 
rises to ~30% at z=l and to ~60% at z=2; for type la SN 
it is about 2% at z=0, 15% at z=l and 35% at z=2. This 
hidden amount of SNe must be taken into account when 
using the SN rate to put constraints on the cosmic SFD 
and to design the future SN searches at high redshift. 
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